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20.  Abstract  (continued) 

^including  both  the  single  event  and  multiple  event  records,  were 
provided  to  Systems,  Science  and  Software  by  the  VELA 
Seismological  Center  (VSC) . 

The  data  were  analyzed  by  the  MARS  (Multiple  Arrival  Recogni- 
tion System)  program  which  uses  a series  of  narrow-band  filters 
to  identify  phase  arrivals.  Their  amplitude  and  arrival  time 
are  accurately  preserved.  The  veuriedjle  frequency  magnitude  (VFM) 
discriminant  is  also  computed  and  the  multiple  explosion  records 
are  clearly  identified  as  being  from  an  explosion. 

For  decomposition  of  the  multiple  event  the  MARS  output  for 
the  single  and  multiple  event  records  were  cross-correlated.  Two 
main  groups  of  events  separated  by  6.7  to  7.8  seconds,  depending 
on  azimuth,  were  identified.  Each  of  these  groups  appeared  to 
include  two  or  more  explosions. 

Subsequent  comparison  with  the  actual  array  configuration 
showed  that  the  two  main  groups  of  events  had  been  correctly 
identified.  There  were  two  events  in  each  group  and  their  lag 
times  were  correctly  identified  within  0.1  second.  For  these 
four  events  the  relative  em^litude  estimates  were  correct  within 
10-20  percent.  The  actual  array  also  included  two  other 
smaller  and  later  events  that  were  not  identified  by  our 
analysis. 


UNCLASSIFIED 

(■CUMITV  CLAturiCATtON  OF  THIS  FAOtdntM  Dmm  tnfnd) 


TABLE  OP  CONTENTS 


INTRODUCTION  AND  SUMMARY 


DATA 


EVENT  IDENTIFICATION  USING  THE  VARIABLE 
FREQUENCY  MAGNITUDE 


DECOMPOSITION  OF  THE  MULTIPLE  EVENT 


4.1  INTRODUCTION 


4.2  IDENTIFICATION  OP  EVENTS 


4.3  AUTOMATION  OP  EVENT  IDENTIFICATION 
BY  CROSS-CORRELATION 


4.4  AMPLITUDE  OF  THE  INDIVIDUAL  EVENTS 


4.5  ORIENTATION  OF  THE  ARRAY 


COMPARISON  TO  ACTUAL  ARRAY  CONFIGURATION 


LIST  OF  ILLUSTRATIONS 


The  data  for  Station  ANNO  are  shovm 


The  data  are  shown  for  Station  MAIO 


The  variable  frequency  magnitude  estimates 
[mb(f)]  from  the  seismograms  of  Figures  1 
and  2 are  compared  to  those  from  a popula- 
tion of  Eurasian  earthquakes  emd  presumed 
explosions  recorded  at  LASA  

The  variable  frequency  magnitude  estimates 
from  the  seismograms  of  Figures  1 and  2 
are  compared  to  those  from  a population  of 
Gulf  of  California  earthquakes  and  NTS 
explosions  recorded  at  YKA  

The  relative  peak  amplitudes  are  plotted 
as  functions  of  frequency  and  arrival  time 
for  Station  ANNO 


The  relative  peak  amplitudes  are  plotted 
as  functions  of  frequency  and  time  for 
Station  MAIO  


The  cross-correlation  function  is  plotted 
for  direct  correlation  of  the  single  and 
multiple  event  records  for  MAIO 

The  cross-correlation  functions  for  ANMO 
are  shown  for  correlation  of  the  narrow- 
band  filtered  signals  from  the  single  and 
multiple  events.  The  correlation  is  done 
filter-by-filter  and  the  results  are  shown 
for  five  center  frequencies  (f^)  

Results  analogous  to  those  in  Figxire  8 are 
shown  for  MAIO  


The  summation  of  the  cross-correlation  func 
tions  of  Figure  8 (AMMO)  is  shown 


The  cross-correlation  function  for  MAIO  is 
shown  in  the  same  form  as  that  in  the 
previous  figure 


12.  The  geometry  is  shown  for  two  stations 

recording  two  events  detonated  in  sequence, 
first  A,  then  b 31 


first  A,  then  b 31 

* 

13.  The  actual  array  configuration  is  shown  ....  33 

14.  The  summation  of  the  cross-correlation 

functions  of  Figure  8 (ANNO)  is  shown 38 

^ 15.  The  cross-correlation  function  for  MAIO  is 

shown  in  the  same  form  as  that  in  the 

previous  figure  39 


9 


9 


LIST  OF  TABLES 

Table  Page 

1.  Individual  Events  In  the  Multiple 

Explosion  Array  3 

2.  Suimnary  of  Event  Lag  Times 18 

3.  Summary  of  Event  Lag  Times 27 

4.  Normalized  An^litudes  of  the  Individual 

Events  in  the  Multiple  Explosion  Array 29 

5.  Data  For  Construction  of  Multiple  Event 

Seisinogreuns 34 

6.  Comparison  of  Actual  and  Computed  Data 

for  ANMO 35 

7.  Comparison  of  Actual  emd  Computed  Data 

for  MAIO 36 


I.  INTRODUCTION  AND  SUMMARY 


f 

I 

L 

I 


f 

« 


» 


f 


t 


I 


L 


The  objective  of  this  study  is  to  detect  amd  identify 
the  individual  events  in  a hypothetical  multiple  explosion. 

The  data  are  simulated  far-field  short  period  body  wave  record- 
ings of  that  event.  Seismograms  were  provided  in  digital  form 
for  two  stations  located  about  19  and  91  degrees  from  the 
event  and  at  azimuths  of  N127"W  emd  N9*E.  The  two  multiple 
event  records  were  presumably  constructed  by  lagging < scaling 
and  summing  a single  event  seismogram  for  the  particular 
source  region-receiver  combination  of  interest.  The  single 
event  record  was  also  provided  in  digital  form. 

The  signal  analysis  program  called  MARS  (Multiple 
Arrival  Recognition  System)  was  used  to  decompose  the  multiple 
event.  In  this  program  the  seismograms  are  filtered  by  a 
series  of  narrow-band  filters  to  identify  phase  arrivals. 

The  amplitude  and  arrival  time  (in  a narrow  frequency  band) 
of  the  identified  phases  are  accurately  preserved. 

A byproduct  of  the  analysis  is  the  variable  frequency 
magnitude  [m^(f)]  for  each  seismogram.  The  ^(f)  values  form 
the  basis  for  a powerful  earthquake/explosion  discriminant 
(Savino  emd  Archambeau,  1974) . The  first  question  addressed 
is,  do  the  seismograms  appear  to  be  from  an  explosion?  This 
question  is  addressed  in  Section  III.  We  find  that  all  four 
seismograms,  the  two  single  and  the  two  multiple  event  records, 
are  clearly  recordings  of  an  explosion  event. 

Decomposition  of  the  multiple  event  record  basically 
requires  a cross-correlation  of  the  pattern  of  phase  arrivals 
in  the  single  event  records  with  those  on  the  multiple  event 
records.  In  Section  IV  we  describe  this  analysis  and  give  the 
results.  The  correlation  is  first  done  by  eye  and  then  is 
done  automatically  by  cross-correlating  the  narrow-band  filter 
output  filter-by- filter  and  summing  the  resulting  correlation 
functions • 
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hint  of  Event  6 and  no  indication  of  the  presence  of  Event  7. 
The  latter,  we  are  fairly  sure,  is  a misidentification  (in 
the  correlation  by  eye)  that  arises  due  to  interfering  phases 
from  earlier  events  in  the  sequence.  The  Event  6 may  or  may 
not  be  present.  Our  best  estimate  is  that  i^  probably  is  not. 

Comparison  to  the  Correct  Answers 

After  completing  the  first  draft  of  the  report  which 
is  summarized  in  the  preceding  paragraphs,  we  were  provided 
with  the  correct  answers.  These  are  compared  to  our  estimates 
in  Section  V.  Briefly  summarized,  the  results  from  the  com- 
parison are  these.  There  were  six  events  with  relative  ampli- 
tudes 0.4,  0.4,  0.8,  0.4,  0.2,  respectively.  Our  automated 
cross-correlation  correctly  identified  the  first  four  of  these 
with  a maximum  error  in  the  lag  time  of  less  than  0.1  seconds. 
For  these  four  our  amplitude  estimates  were  correct  within  10- 
20  percent.  The  last  two  events  were  entirely  missed,  pre- 
sumably because  they  are  smaller. 

We  also  ventured  an  estimate  of  the  array  orientation 
which  can  be  compared  to  the  actual  orientation  shown  in 
Figure  13.  As  was  thought,  the  first  two  events  proceeded 
(generally)  from  south-to-north.  Also,  the  first  event  in 
the  second  group  is  12  km  south  of  the  last  event  in  the 
first  group,  which  is  close  to  our  estimate. 

The  best  results  were  achieved  with  the  automated 
cross-correlation,  which  is  encouraging.  The  more  arduous 
Informal  correlation  misidentlfied  several  events  and  gave 
relatively  poor  estimates  of  the  relative  amplitudes.  How- 
ever, this  method  was  helpful  in  guiding  our  thinking  and 
in  confirming  the  more  abstract  results  from  the  automated 
cross-correlation . 


The  data  used  for  this  test  consist  of  short-period 
vertical  component  recordings  from  stations  ANMO  emd  MAIO  for 
two  events.  The  first  is  a single  explosion  at  the  Soviet 
test  site  near  Semipalatinsk.  The  second  event  is  a synthetic 
multiple  explosion  constructed  by  scaling,  lagging  and  summing 
the  single  event  records.  The  four  records  (two  for  ANMO  and 
two  for  MAIO)  are  shown  in  Figures  1 and  2.  All  data  were  pro' 
vided  by  the  SDAC  in  digital  form  with  a sampling  rate  of  20 
points  per  second.  The  station  ANMO  is  about  91  degrees  from 
Semipalatinsk  at  an  azimuth  of  9 degrees  east  from  north. 
Statio^n  MAIO  is  at  an  epicentral  distance  of  about  19  degrees 
and  an  azimuth  of  233  degrees  east  from  north.  The  two  sta- 
tions are  therefore  136  degrees  apart. 


5 


) 

I 

5 


f 

I III.  EVENT  IDENTIFICATION  USING  THE  VARIABLE 

' ’ FREQUENCY  MAGNITUDE 

An  obvious  question  one  might  ask  about  the  simulated 
multiple  explosion  under  study  is,  cam  the  multiple  explosion 
be  clearly  Identified  as  an  explosion  event  rather  than  an 
earthquake?  Since  the  n^(f)  values  used  for  the  VFM  (Vari- 
able Frequency  Magnitude)  discrimln2uit  (Savlno  and  Archambeau, 
1974)  are  byproducts  of  our  narrow-band  filter  analysis,  we 
can  plot  these  values  together  with  earthquaUce  and  explosion 
data  previously  analyzed. 

In  Figure  3 is  shown  a plot  of  m^  (0.45  Hz)  versus 
(2.25  Hz)  for  LASA  recordings  of  shallow  Eurasian  earthquakes 
and  presumed  explosions.  Also  shown  are  the  same  quantities 
for  several  simulated  multiple  explosion  seismograms  that  were 
constructed  in  much  the  same  way  as  those  being  studied  here 
(for  a detailed  description  of  Figure  3 see  Chapter  VI  of 
Bache,  ^ al . , 1976).  These  earlier  multiple  explosion 
scenarios  were  specifically  designed  to  appear  earthquake- 
like in  terms  of  the  complexity  and  first  motion  dis- 

criminants. These  simulated  multiple  explosions  clearly 
separate  from  the  shallow  earthquakes  to  which  they  are  com- 
pared in  Figure  3. 

The  station  ANMO  is  at  a distance  (s91  degrees)  from 
the  Semipalatlnsk  event  that  is  about  the  same  as  most  of  the 
events  in  the  LASA  data  plotted  in  Figxire  3.  We  plot  the  data 
from  both  the  single  event  (Z)  and  the  simulated  multiple 
event  (ZSUM)  on  Figure  3 and  see  that  they  are  both  clearly 
I identified  as  explosions.  We  point  out  that  we  are  uncertain 

I of  the  static  Instrument  gain  for  the  records  and  it  is  pos- 

sible that  a scale  factor  should  be  added  to  all  ^(f)  values, 
I ' moving  them  along  a 45  degree  line. 


45  Hz) 


^ (2.25  Rz) 


Figure  3.  The  variable  frequency  magnitude  estimates  [^(f)] 
from  the  seismograms  of  Figures  1 and  2 are 
conqpared  to  those  from  a population  of  Eurasian 
earthquakes  and  presumed  explosions  recorded  at 
LASA. 


The  ^(f)  values  for  station  MAIO  are  also  plotted  on 
Figure  3 and  fall  well  outside  all  other  values.  This  is  pre- 
sumably because  MAIO  is  at  a much  closer  distance  (»19  degrees) . 

For  a more  reasoneUdle  comparison  to  the  n^(f)  values 
for  MAIO,  we  select  a previously  analyzed  population  of  North 
American  events  recorded  at  the  Cemadian  array  YXA.  These 
data  are  shown  in  Figure  4,  together  with  the  ^(f)  values 
from  MAIO  and  ANMO.  Since  we  do  not  know  the  static  instru- 
ment gains,  these  values  can  move  along  a 45  degree  line  as 
indicated  in  the  figure.  Once  again,  both  the  single  and 
multiple  explosion  events  are  clearly  identified  as  explosions. 


Figure  4.  The  varieUble  frequency  magnitude  estimates  from 
the  seismograms  of  Figures  1 and  2 are  compared 
to  those  from  a population  of  Gulf  of  California 
earthquakes  and  NTS  e3q>losions  recorded  at  YKh. 


S 


11 


i 

I 

I 


1 


I 


IV.  DECOMPOSITION  OF  THE  MULTIPLE  EVENT 

4.1  INTRODUCTION 

The  method  used  for  identifying  the  individual  events 
in  a multiple  explosion  array  was  described  in  some  detail  in 
an  earlier  report  (Leutibert,  ^ al . , 1977).  In  the  earlier 
work  our  attention  was  directed  almost  entirely  to  near-field 
data.  In  the  current  study  we  are  concerned  with  recordings 
at  large  ranges  (approximately  19  and  91  degrees)  where  the 
high  frequency  information  has  been  greatly  attenuated. 

Briefly  described,  the  method  ezaployed  is  as  follows. 

Each  seismogram  to  be  analyzed  is  Fourier  trzmsformed  and 
filtered  by  a narrow-band,  Gaussian  filter.  The  Hilbert  trams- 
form  of  the  filtered  spectrixm  is  then  constructed  and  inverse- 
transformed  to  obtain  the  envelope  function.  The  peaks  of  the 
envelope  function  Indicate  phase  arrivals  of  energy  in  a 
narrow  b^uld  of  frequencies  near  the  filter  center  frequency. 

The  arrival  times  are  accurately  preserved  by  the  times  of 
these  peaks.  Further,  the  relative  amplitudes  of  the  peeUcs 
reflect  the  relative  amplitudes  of  the  arriving  phases. 

For  the  experiment  vdiich  this  report  describes  we  have 
a multiple  event  record  from  two  stations.  These  multiple 
event  records  are  presumed  to  have  been  constructed  by  lagging, 
scaling  and  summing  a single  event  record  according  to  some 
formula  based  on  the  timing  and  spacing  of  the  individual  | 

events  in  the  hypothetical  multiple  explosion  array.  Me  also  | 

have  been  provided  with  the  two  single  event  records. 

In  order  to  decompose  the  multiple  event  record  we  do  j 

the  following.  All  four  records  are  filtered  by  narrow-band 
filters  at  series  of  center  frequencies  spanning  the  range  at 
^Ich  spectral  energy  is  present.  Our  task  is  to  Identify  the 
arriving  phases  on  the  single  event  records  and  then  to 
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correlate  this  pattern  with  the  identified  phase  arrivals  on 
the  multiple  event  records  to  Identify  the  repeated  presence 
of  the  single  event. 

What  are  the  constraints  on  this  method?  In  our  pre- 
vious study  (Lambert,  et  , 1977),  we  found  that  we  were 
able  to  separate  individual  arrivals  whenever  there  was 
"sufficient”  signal  energy  at  frequencies  greater  them  about 
3.5  times  the  inverse  of  the  lag  time  between  these  arriving 
signals.  For  teleseismic  short  period  recordings  this  means 
that  identifying  events  separated  by  less  than  a second  can 
only  be  achieved  if  there  is  considerable  signal  energy  at 
frequencies  of  3.5  Hz  and  higher.  Generally  this  is  not  the 
case. 

Basically  our  problem  is  one  of  cross-correlation  of 
a particular  kind  of  output.  In  Section  4.2  we  describe  the 
results  achieved  by  carrying  out  the  cross-correlation  by 
eye.  While  it  may  work  best  because  of  the  opportunity  to 
inject  human  judgement,  this  is  certainly  a slow  emd  arduous 
way  to  proceed.  Therefore,  the  cross-correlation  was  auto- 
mated in  a manner  described  in  Section  4.3.  The  automated 
method  worked  very  well  emd  gives  the  best  results.  The 
informal  correlation  is  described  because  of  the  insight  it 
gives  to  the  procedure. 

4.2  IDENTIFICATION  OF  EVENTS 

A series  of  twenty  narrow-band  filters  with  center 
frequencies  ranging  from  0.3  to  9.0  Hz  were  applied  to  the 
single  and  multiple  event  signals  (Figures  1 emd  2)  from 
the  stations  MAIO  and  ANMO.  A parameter  of  some  minor  impor- 
tance is  the  filter  Q or  width.  After  several  experiments, 
a value  that  seemed  to  be  optimum  was  selected.  The  results 
were  then  plotted  as  shown  in  Figures  5 and  €.  For  each  seis- 
mogram these  are  plots  of  the  time  of  arrival  of  the  largest 
envelope  peaks  versus  the  filter  center  frequency  (f^) . 

13 
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Figure  5.  The  relative  peak  an^litudes  are  plotted  as  functions  of  frequency  and  arrival 
tisie  for  Station  ANMO.  The  plot  for  the  single  event  record  is  on  top  with 
that  for  the  multiple  event  below. 
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For  each  the  plot  symbols  are  numerals  that  indicate  the 
relative  amplitudes  of  the  peaks  compared  to  the  largest  for 
that  frequency,  which  is  indicated  by  an  asterisk.  A "9" 
indicates  that  the  amplitude  is  greater  than  90  percent  of 
the  largest,  etc.,  down  to  "0”  which  indicates  peaks  with 
amplitudes  from  1 percent  to  10  percent  of  the  largest. 

The  next  step  is  to  identify  the  nondispersed  phase 
arrivals.  These  are  arrivals  for  which  there  are  envelope 
peaks  of  roughly  the  same  relative  size  persisting  over  a 
bemd  of  frequencies.  We  began  by  identifying  these  by  eye, 
a task  that  requires  some  judgement,  but  one  that  is  neces- 
S2u:y  before  attempting  to  automate  the  process.  The  results 
obtained  by  direct  examination  of  the  plots  in  Figures  5 euid 
6 are  described  in  this  section  while  a later  Section  (4.3) 
is  devoted  to  the  results  from  an  automated  analysis. 

In  Figures  5 and  6 solid  lines  are  drawn  through  the 
major  phase  arrivals.  The  dashed  lines  indicate  distinct 
arrivals  that  may  be  present  but  are  less  certainly  identified. 

To  decompose  the  multiple  event  records  the  task  is 
basically  to  cross->correlate  the  two  pairs  of  single  and 
multiple  event  records  in  Figures  5 and  6.  For  each  we  then 
have  "single  station"  decompositions.  The  results  from  the 
two  can  then  be  compared  to  in^rove  the  resolution.  However, 
a complication  is  that  stations  at  different  azimuths  will 
generally  not  see  the  seune  lag  times  (recall  from  Section  II 
that  ANMO  and  MAIO  eure  136  degrees  apart  in  azimuth) . 

For  Station  ANMO  a comparison  of  the  single  and  multiple 
event  plots  clearly  indicates  the  presence  of  two  major  groups 
of  events  separated  by  8 seconds.  There  are  indications  of 
at  least  two  explosions  in  each  group  with  the  four  identified 
events  occurring  at  0.0,  0.81,  7.80  and  8.81  seconds.  There 
is  a less  certain  Indication  of  a later  event  separated  by 


14.67  seconds  from  the  first,  but  this  may  be  due  to  construc- 
tive interference  of  later  arrivals  from  the  first  four  events. 

The  Station  MAIO  is  much  closer  to  the  source  than 
ANMO  (19  degrees  versus  91  degrees)  and  therefore  contains 
considerably  more  signal  energy  at  high  frequencies.  This 
is  apparent  in  the  envelope  peak  plots  for  MAIO  in  Figure  6. 
Once  again,  the  single  and  multiple  event  records  are  com- 
pared to  discover  repeated  patterns.  Individual  events  can 
tentatively  be  identified  at  0.0,  0.51,  0.9,  6.7,  7.6,  13.5 
and  21.3  seconds.  The  events  at  0.51  and  7.6  seconds  are 
relatively  less  clear  them  those  at  0.0,  0.9  and  6.7  seconds. 
Those  at  13.5  and  21.3  seconds  may  not  be  real  but  may  appear 
due  to  interference  from  later  arrivals  associated  with  the 
earlier  events. 

To  obtain  additional  confidence  in  our  event  decon^si- 
tion,  we  compare  results  from  the  two  stations.  These  are 
summarized  in  Table  2.  In  order  to  combine  the  information 
from  the  two  stations,  we  must  make  some  assumptions  about  the 
origin  of  the  lag  time  between  events.  There  are  basically 
two  ways  for  these  lag  times  to  arise.  First,  the  explosions 
could  be  closely  spaced,  but  detonated  in  some  sequence. 

Second,  the  explosions  could  be  detonated  slmult2meously  and 
the  lag  time  be  due  to  them  being  varying  distemces  from  the 
receiver.  However,  at  91  degrees  the  dT/dA  is  about  4.7 
seconds/degree  and  lag  times  of  a second  require  differential 
path  lengths  of  about  23  km.  At  19  degrees  the  event  separa- 
tion per  second  of  time  lag  is  about  9 km.  Therefore,  it  is 
likely  that  most  of  the  time  lag  between  events  (certainly 
the  larger  lag  times)  is  due  to  them  being  detonated  in 
sequence . 

Comparing  the  event  arrival  times  in  Table  2,  we  see 
three  possible  groups  of  events.  Events  1,  2 and  3 at  MAIO 
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TABLE  2 


SUMMARY  OF  EVENT  LAG  TIMES 


Identified  Event 


ANMO 

Lag  Time  (sec) 


MAIO 

Lag  Time  (sec) 


0.0* 


0.0* 


0.51* 


0.81 


0.90 


7.80 


6.70 


8.81* 


7.60* 


14.67** 


13.51** 


21.31**  - 


* Less  correlation  than  those  without  asterisks. 

**  These  correlations  could  be  due  to  constructive 
Interference  from  earlier  events. 


I 


i 


and  1 and  3 at  ANMC  comprise  the  first  group.  Lag  time  varia*  | 

tlons  between  stations  of  th^  order  represented  could  easily  be 
due  to  azimuthal  effects.  The  second  group  of  events  are  those 
indicated  as  Events  4 emd  5.  There  is  then  a single  event 
identified  as  6 at  both  stations.  Finally,  there  is  the  event 
at  21.31  on  the  MAIO  record  for  which  no  corresponding  arrival 
is  seen  at  ANMO.  This  is  probably  a spurious  event  identifi- 
cation. 

\ 

In  summary,  two  groups  of  explosion  events  are  clearly 
present,  being  separated  in  time  by  about  6 to  9 seconds. 

Actually,  this  is  fairly  clear  from  direct  examination  of  the 
original  seismogrcuns  (Figures  1 and  2) . Each  of  the  two  main 
groups  seems  to  be  composed  of  two  or  more  explosions  with 
separation  times  on  the  order  of  a second.  Finally,  there 
may  be  at  least  one  other  event  lagged  some  5 to  6 seconds 
after  the  second  group. 

Small  differences  in  lag  times  between  events  can 
potentially  be  used  to  ascertain  the  orientation  of  the 
multiple  explosion  array  if  we  assume  it  to  be  linear.  This 
subject  is  discussed  in  Section  4.5. 

4.3  AUTOMATION  OF  EVENT  IDENTIFICATION  BY  CROSS -CORRELATION 

The  event  identification  results  of  the  previous  section 
were  obtained  by  a careful  ex2unination  of  the  narrow-band  filter 
output.  While  slow  and  arduous,  this  method  gives  maximum  rein 
to  human  judgement  and  intuition.  As  was  pointed  out,  the 
event  identification  process  is  essentially  one  of  cross-  I 

correlation  using  the  single-  and  multiple-  event  records.  In  | 

this  section  we  describe  the  results  obtained  when  this  process 
is  automated. 

The  cross-correlation  of  two  signals  is  defined  by 
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c(t)  - ^ fJ((d)  d w, 
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where  F2(u)  is  t^e  complex  conjugate  of  the  Fourier  transform 
of  F2 (t) . The  cross-correlation  results  presented  in  this 
report  are  from  an  Systems,  Science  and  Software  (S*)  modified 
version  of  the  program  COLLAPSE  (Alexander  and  Lambert,  1971) . 

To  decompose  the  multiple  event,  one  of  the  simplest 
things  one  could  try  would  be  to  cross-correlate  the  time  sig- 

I nals  for  the  single  and  multiple  events  at  each  station.  The 

results  for  Station  MAIO  are  shown  in  Figure  7.  Only  the  two 
ma!in  groups  of  events  separated  by  about  six  seconds  can  be 

I seen  and  even  these  are  not  indicated  unambiguously.  Elemen- 

tary cross-correlation  of  the  signals  clearly  gives  insufficient 
resolution. 

In  view  of  our  analysis  of  the  data  in  Figures  5 and 
6,  a better  way  to  proceed  is  to  cross-correlate  the  narrow- 
band  filter  output.  This  is  done  filter-by-filter  and  the 
results  are  summed  to  give  the  final  cross-correlation  func- 
tion. In  Figures  8 and  9 we  show  the  cross-correlation  func- 
tions for  each  station  for  a selected  set  of  center  frequencies. 
For  the  individual  filters  only  the  two  main  groups  of  events 
are  clearly  identified  at  about  14  and  21  seconds  into  time 
signal.  At  this  stage  the  narrow-band  filter  output  cross- 
correlation gives  results  little  better  than  those  from  direct 
cross-correlation  of  the  seismograms  (Figure  7) . However,  the 
results  are  considerably  enhanced  by  summing  the  individual 
filter  correlations. 

In  Figures  10  euid  11  we  show  the  final  cross-correlation 
functions  for  ANMO  and  MAIO.  This  function  is  a sum  of  the 

f 

I correlation  functions  for  the  individual  filters  (Figures  8 and 
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Phase  arrival  times 
from  Table  2 


1.0000*02 


Figure  7.  The  cross-correlation  function  is  plotted  for  direct 
correlation  of  the  single  and  multiple  event  records 
for  NATO. 


Figure  8.  The  ANMO  cross-correlation  functions  are  shown  for 
correlation  of  the  narrow-band  filtered  signals 


from  the  single  and  multiple  events.  The  correla- 
tion is  done  filter-by-fllter  and  the  results  are 
shown  for  five  center  frequencies  (f  ) . The  phase 
arrival  times  from  Table  2 are  Indicated  on  each 
plot. 
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sun  OF  CONVOLUTIONS  OF  FILTER  18^012  UITH  SIGNAL  18^012 


0.0000  S.OOOO'OO  1.C>000*01  1.S00O«01  2.0000*01  2.5000*01  3.0000*01  3.5000*01  4.0000*01 


TinE  (SEC) 


Figure  10.  The  summation  of  the  cross-correlation  functions  of 
Figure  8 (ANMO)  Is  shown.  The  four  Identifiable 
arrivals  listed  In  Table  3 are  Indicated  by  arrows. 
Also  Indicated  are  the  arrival  times  from  Table  2 
that  do  not  correspond  to  Identifiable  peaks  In  the 
cross-correlation  function. 


sun  OF  CONVOLUTIONS  OF  FILTER  18^032  UITH  SIGNAL  18^032 


Figur*  II.  The  cross-correlation  function  for  MAIO  is  shown  in 
the  sane  form  as  that  in  the  previous  figure. 


9).  The  answer  depends  to  some  extent  on  the  selection  of 
filters  making  up  the  sum.  The  particular  combination  shown 
seemed  to  be  the  best  of  those  tried. 

The  individual  event  arrivals  identified  by  the  cross- 
correlation functions  are  indicated  in  Figures  10  and  11  and 
are  tabulated  in  Table  3.  Also  summarized  are  the  arrival 
times  from  the  analysis  of  Section  4.2.  For  the  four  events 
in  common,  the  cross-correlation  results  seem  to  agree  rather 
well  with  those  from  direct  examination  of  the  narrow-band 
filter  output.  Taken  together,  the  two  give  a more  complete 
picture  of  the  event. 

We  conclude  that  there  are  two  main  groups  of  events 
(1-3  and  4-5)  separated  by  6.7  to  7.8  seconds,  depending  on 
the  azimuth.  Each  group  consists  of  two  or  more  events 
separated  by  lag  times  small  enough  to  be  at  the  limits  of 
our  resolving  power.  The  cross-correlation  fxinctions  seems 
to  indicate  that  Event  7 is  almost  certainly  not  present 
2uid  was  misidentified  by  the  analysis  of  Section  4.2.  The 
Event  6 is  difficult  to  judge.  There  is  a hint  of  something 
in  the  cross-correlation  function,  but  it  is  certainly  very 
weak.  Our  best  judgement  is  that  this  event  is  not  present. 


4.4  AMPLITUDE  OF  THE  INDIVIDUAL  EVENTS 

As  with  the  lag  time  determination,  we  now  estimate 
the  amplitude  of  the  individual  events  with  two  different 
methods.  We  begin  by  using  the  data  plotted  in  Figures  5 
and  6.  Each  event  arrival  is  identified  by  envelope  peaks 
that  arrive  at  the  same  time  over  a band  of  frequencies. 

The  amplitude  of  each  envelope  peak  accurately  reflects  the 
energy  of  the  phase  arrival  at  this  frequency.  A particular 
event  arrival,  say  that  at  27.5  seconds  on  the  MAIO  plot  of 
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(azimuth  N127'*W) . Since  the  dA/dT  is  about  23  km/sec  at  ANMO 
versus  9 km/sec  at  MAIO,  a linear  array  in  line  with  ANMO 
requires  the  smallest  shot  separation  consistent  with  these 
lag  times. 

If  any  two  events  were  known  to  have  occurred  simul- 
temeously  and  the  differential  lag  times  were  known  exactly, 
we  could  fix  the  array  orientation.  However,  the  lag  times 
we  have  determined  are  uncertain  by  a few  tenths  of  a second 
(see  the  comparison  in  Table  3)  and  it  is  differences  of  this 
order  that  are  indicative  of  the  orientation.  Further,  since 
the  events  may  occur  in  sequence,  the  problem  is  underdeter- 
mined with  two  stations.  That  is,  for  any  pair  of  stations 
recording  a pair  of  events  we  can  write  the  set  of  equations 

ti  » Pj^D  cos  <l>i  + Tq  / 

^2  * ^2°  *^2  ^o  ' 

♦2  ” ^1  “ • 

where  t^  and  t^  are  the  lag  times  perceived  at  Stations  1 and  2, 
p^  and  P2  are  the  ray  parameters  for  these  stations,  T^  is  the 
detonation  lag  time  and  the  other  quantities  are  defined  in 
Figure  12.  These  are  three  equations  in  the  four  unknowns  T^, 

D,  emd  ^2*  events  were  detonated  simultaneously 

(Tq  ■ 0) , these  can  be  solved.  For  T^  7^  0 a third  station  is 
sufficient  to  uniquely  determine  the  solution.  This,  of 
course,  assumes  no  errors  in  the  lag  times.  If  errors  are 
present,  more  stations  are  needed  so  a best  fitting  solution 
to  an  overdeterroined  set  of  equations  can  be  constructed. 

If  the  events  are  in  a linear  array,  it  doesn't  matter 
how  many  there  are.  Each  additional  event  adds  two  equations 
like  the  first  two  listed  above  and  two  unknowns,  the  detona- 
tion lag  time  and  the  separation.  Therefore,  three  stations 
are,  in  principle,  enough  to  uniquely  determine  the  orienta- 
tion of  the  array. 
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V.  COMPARISON  TO  ACTUAL  ARRAY  CONFIGURATION 


* 

I 

After  completing  the  analysis  described  in  Section  IV,  j 

the  details  of  the  actual  arr^y  configuration  were  provided 
for  comparison  to  our  answers  (Captain  M.  J.  Shore,  personal 
communication) . There  were  six  events  in  a linear  array 
oriented  N45*W.  The  detonation  sequence  and  event  spacing 
are  shown  schematically  in  Figure  13. 

The  multiple  event  records  were  constructed  by  lagging, 
scaling  and  summing  individual  event  records  at  the  two 
stations.  It  is  the  time  lag  and  relative  size  of  events  on 
each  record  that  we  determined  in  our  analysis.  The  actual 
numbers  used  to  construct  the  seismograms  are  summarized  in 
Table  5. 

Before  comparing  the  actual  values  to  those  deduced 
in  Section  IV,  we  make  some  comments  about  the  time  lags  in 
Table  5.  The  equations  relating  the  time  lags  to  the  array  i 

configuration  are  given  in  Section  4.5.  If  we  use  the  station 
locations  (distance,  azimuth  with  respect  to  the  event)  given  . 

in  tt^iis  report,  the  time  lags  in  TeUsle  5 are  inconsistent  with  ' 

the  array  configuration,  though  the  deviation  is  small.  There-  ! 

fore,  the  distances  and  azimuths  used  for  computing  the  actual  ' 

lags  must  have  been  different  than  those  mentioned  in  this 
report . 

The  actual  and  computed  lag  times  and  relative  amplitudes 
are  compared  in  Tables  6 and  7.  As  might  have  been  expected, 
the  lag  times  from  the  automated  cross-correlation  of  Section 
4.3  are  the  most  accurate.  With  this  method  we  correctly 
identified  Events  1,  2,  3 and  4 at  both  stations  with  the 
lag  time  error  being  less  than  0.1  seconds  in  every  instance. 

The  Events  5 and  6 are  relatively  small  and  were  not  Identified 

by  our  analysis.  I| 
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TABLE  5 

DATA  FOR  CONSTRUCTION  OF  MULTIPLE  EVENT  SEISMOGRAMS 
(from  Capt.  M.  J.  Shore^  Personal  Communication) 


Event 

Spacing 
from  1 (km) 

Relative 

Amplitude 

ANNO 

Time  Lag  (sec) 

MAIO 

Time  Lag  (sec) 

1 

0.0 

0.4 

0.0 

0.0 

2 

4.0 

0.4 

0.38 

0.94 

3 

-8.0 

0.8 

7.94 

6.82 

4 

-6.5 

0.4 

8.69 

7.79 

5 

0.25 

0.2 

10.49 

10.53 

6 

6.0 

0.2 

10.32 

11.16 

1 
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Figure  14. 


The  summation  of  the  cross-correlation  functions  of 
Figure  8 (ANMO)  is.  shoim.  The  four  Identifiable 
arrivals  listed  in  Table  3 are  indicated  by  arrows. 

We  Indicated  the  arrival  times  from  Table  2 that  do 
not  correspond  to  Identifiable  peaks  in  the  cross- 
correlation  function  and  also  those  that  arm  actually 
present  (Personal  Communication,  1977)  but  not 
identified; 
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